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While some experiments seem to confirm the contention that the axial 
buckling load of a pressurized cylinder approaches the classical value (i.e., 
the value predicted by the linear buckling equations) for sufficiently large 
internal pressure, other tests indicate a much smaller buckling load increase 
resulting from pressurization. Here it is shown that the buckling load of 
an elastic shell with asymmetric imperfections|but sufficiently free of 
axisymmetric imperfections closely coincides with the classical value for 

relatively small values q f internal pressure; while the ^buckling load of a 

y 

shell with a predominance of axisymmetric imperfection can remain well below 
the classical value for the entire range of internal pressures for which the 
shell buckles elastically. Of particular interest is the/ calculation of an 
upper bound to the buckling load as predicted by the non-linear Donnell-shell 
equations for a shell with axisymmetric imperfections.;' 



This work was supported by the National Aeronautics and Space Administration 
under Contract No. NsG-559 with Harvard University. 

GPO PRICE $ 

OTS PRICE(S) $ 


Hard copy (HC) 
Microfiche (MF) 


s£jv 
■ Sd 


2 


INTRODUCTION 

The several published analyses of pressurized cylindrical shells under 
axial compression do not adequately explain the variety of behavior which has 
been reported for such structures. The linear buckling equations predict that 
the buckling parameter 

fjCc-f) 

is unity for pressurized and unpressurized shells. Here o is the applied 
compressive stress and pR/2h is the axial stress resulting from pressurization 
as depicted in Figure 1. Experiments indicate that this buckling parameter is 
usually on the order of one half or one third for unpressurized shells and is 
larger for pressurized shells. In some tests the parameter is unity for shells 
under sufficient internal pressure, while in other cases this parameter remains 
well below unity for the entire range of pressures for which the shell buckles 
elastically. 

j_iu, ul due duu ouiwdLL^ v . ) ubeu non-unear Ducxung equations ror a pertect 

shell and employed the energy criteria of buckling to show that the buckling 

parameter, as defined here, increases from .62 at p = 0 to unity at 
2 2 

pR /Eh = .17 . Thieleman (2) also reported calculations for initially perfect 
shells based on a non-linear analysis. He obtains load-deflection curves from 
which is obtained the minimum load the shell can support following buckling. 

This minimum load increases with increased internal pressure. 

The role of shell imperfections, known to be the main degrading factor in 
unpressurized shells, was not considered in either of the above mentioned papers. 
Lu and Nash (3) have studied the effect of initial imperfections on the minimum 
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load the shell can support in the post-buckling region. This analysis also 
shows a larger minimum load for pressurized shells than unpressurized ones. 

From a design stand-point, the maximum load the shell can support prior to 
buckling is of more interest than the minimum post-buckling load. Indeed, it is 
this value which most experimenters record and which we will designated as the 
buckling load in this report. Here, approximate solutions to the non-linear, 
Donnell-type shell equations will be obtained which display the role of 
axisymmetric and asymmetric imperfections on the buckling load of a pressurized 
elastic shell. It will be seen that axisymmetric imperfections of a certain 
wavelength are particularly degrading. An upperbound to the exact buckling load 
as predicted by the non-linear equations is obtained for the case of axisymmetric 
imperfections. Finally, these solutions are compared with some experimental 
results, and it is noted that they seem to account for the variety of reported 
behavior. 


DONNELL SHELL EQUATIONS 

xut; eiasuL sneii j_t> maidLteLiieu uy jlls> Lciuj.cti uibpicicemeiit w ^pubiuve 

outward) and an Airy stress function F which gives the resultant membrane 

stresses as N = F . N = F and N = -F . The shell is assumed 
xx ,yy 7 yy ,xx xy ,xy 

to have an initial imperfection in the form of an initial, stress-free radial 
displacement W q . With x and y as the cartesian coordinates in the axial 
and circumferential directions, the equations are 
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where 



is the two dimensional biharmonic operator. 
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CLASSICAL BUCKLING 

We consider a perfect cylinder (W q = 0) under internal pressure p and 
axial compressive stress o (in addition to the axial stress resulting from 
the internal pressure). The state of deformation can be written as 


W | (oh + ^) + pR + w 
F = ( - j oh + \ pR)y 2 + \ pRx 2 + f 


(3) 

( 4 ) 


where the terms added to w and f constitute the prebuckling solution for the 
perfect shell. The classical buckling equations are obtained by substituting 
(3) and (4) in the Donnell-shell equations and then linearizing the resulting 
equations with respect to w and f . The linear buckling equations are 
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Solutions to these equations are well known. * The eigen-value 


x = !Vi = lf (<* 2 + P 2 ) 2 + i— 
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corresponds to a radial deflection mode of the form 


w = cos(Ctq o x/R)cos(Pq o y/R) 


( 8 ) 


and the associated stress function 


f = - 


a 


ERh 

2 . 2 , „ 2. 2 

q o (O + P ) 


cos (aq^x/R) cos (Pq^y/R) 
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where 


4 - 2 '" 2 '’ 2 and p = ^3(1 - v 2 ) 


q^ = 12(1 - V )R /h 


Eh 


The classical buckling parameter for the unpressurized shell, obtained by 
minimizing X as given by equation (7) with respect to a and P , is 

2 

( 10 ) 
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In this paper we will assume the shell is sufficiently long to justify 
neglecting the boundary conditions at the ends of the shell. For certain 
end conditions Hoff (4) has produced a solution to the linear buckling 
equations which predicts buckling at X = 1/2 for an unpressurized 
cylindrical shell. It is not known if the post-buckling behavior in this 
case is unstable, nor is it clear that the shell is significantly imperfection 
sensitive for buckling in this mode. Although in certain instances the shell 
buckling may be characterized by the solution obtained by Hoff, we will assume 
that shell imperfections result in the main reduction of the buckling load 
below what is generally accepted as the "classical 11 value, X = 1 . 
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with an infinite number of associated buckling modes such that a and P 
satisfy 

a 2 + p 2 - a = o (li) 

Included in this set of critical buckling modes are the axi symmetric mode 
w = cos(q Q x/R) and the asymmetric mode w = cos( — q o x/R)cos( q^y/R) . 
The classical buckling load (i.e., value of buckling parameter) for the 
pressurized shell is also unity 



however only the axi symmetric mode, w = cosCq^x/R) , is associated with this 
critical value. All other modes of the form of equation (8) are associated with 
eigenvalues larger than unity. Thus, for example, corresponding to 
w = cos( ■- q o x/R)cos( q^y/R) is the eigen-value \ = 1 4* p . 

NON-LINEAR BUCKLING EQUATIONS 

We must turn to the non-linear equations and the effect of initial imper- 
fections to explain the discrepancy between the predictions of the linear or 
classical theory and common experimental observations. It is of particular 
interest that Koiter (5) has presented an analytic procedure for obtaining the 
role of imperfections in imperfection-sensitive structures. This work is perhaps 
more readily available in reference (6) or in a less general form in reference 
(7). His general theory for cylindrical shells (8) indicates that an imperfection 
in the form of the radial displacement component of the axisymmetric buckling 
mode, W =^ccs (q^x/R) , is the most degrading and is able to reduce the buckling 
load of an unpressurized shell under axial compression to one half or even one 
third the classical value for values of M- only a fraction of the shell thickness. 
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In general any radial imperfection pattern of the shell can be represented 
by a double Fourier series in the axial and circumferential coordinates. We 
will restrict ourselves to a consideration of just two terms of such a series - 
one axi symmetric and one asymmetric. Both are taken to be in the form of a 
radial displacement component of a linear buckling mode of the unpressurized 
shell. 

Thus, we study the behavior of a pressurized cylindrical shell under an 
applied axial load where the initial imperfection is assumed to be 

W o = ” h h G 0 S ( < l 0 x / R ) + 1 2 h cos ( \ q Q x/R)cos( J q Q y/R) (12) 

where and f 2 are the ratios of the amplitude of the imperfection to the 

shell thickness. 

Any equilibrium state of the axially loaded cylinder can be written in the 
form of (3) and (4) ; we approximate w by 

11 11 
w = ^h cos(q o x/R) + £ 2 h cos( j q o x/R)cos( j q Q y/R) + l^h sin( j q Q x/R)cos( j q Q y/R) 

(13) 

Here , | 2 an ^ I 3 are th e ratios of the amplitude of the deflection in the 

axisymmetric or asymmetric modes to the shell thickness. 

Solutions to the full non-linear equations (1) and (2) are obtained in the 

% 

following manner. First, since equation (1) is a compatibility^ it is solved 
exactly for F in terms of the assumed W . This is accomplished with the aid 


of equations (3) and (4). Since we are ignoring the end conditions and as the 
average stresses in the shell are given by the polynomial terms in (4), f is 
required to be the sum of terms periodic in the axial and circumferential 
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directions. Secondly, we solve equation (2) (an equilibrium equation) approximately 
by substituting in F and W and then applying the Galkerin procedure. The steps 
of this calculation are easily carried out and the resulting equations are 


1,(1 - X ) + 


h ■ 4 ) 


+ 5 + x h + lk ° 2 


+ i L.#, 

ZOO ^1 


(I 2 + i 2 )(I 2 Si + dj. - 

- i x >s| =0 (14) 


£ 2 (1 + P - X) + | c| 1 S 2 - c^ + cij 2 + (p - \)l 2 + c 2 (| 1 - I 1 )(| 1 I 2 + - i x )) 

+ ^ c 2 (i 2 + i 2 )(g 2 i 2 + \ £ 2 ) - 0 

(15) 

and 

13(1 + P - X) 2 c ^i^3 + c Ml + 25 c ^1 " ^3 + 16 c ^3 = 0 

(16) 

I ^ 

wn ere u - v jvi * v ) • 

A solution to these equations would provide the equilibrium configuration 

of the shell as a function of X, i*e. and as a function of x . 

If X attains a maximum as the compressive axial load is applied then this is 

the value of X associated with the buckling load. It will be denoted by X^ . 

The essential character of the shell behavior in the prebuckiing and immediate 

-2-2 3 

post-buckling description is retained if the terms of order gg , gg and g 
are omitted from the previous three equations. With this simplification the 
equations to be studied are 

£^(1 - X) + -32 c (&2 ~ £ 3 ) "g c ^ 2^2 ~ 


(17) 
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i 2 (l + p - X) +| cg 1 l 2 - cg^ + c g x I 2 - (X - p)I 2 (18) 

and 

S 3 (l + P - X) - f cg^ + = 0 (19) 

For sufficiently small and ^ > these equations yield a buckling 

load which is asymptotic to that predicted by (14), (15) and (16) if p = 0 . 

If p > 0 we cannot make this assertion; however, the upperbound solution for 
the case of axisymmetric imperfections obtained in a later section indicates that 
these equations provide a sufficiently accurate estimate of the buckling load for 
the purposes of this paper. 

The behavior of the perfect shell as predicted by these equations is shown 
in Figure 2. A perfect shell suffers no deformation in the buckling modes 
(ll^Z^s^O) until X reaches unity. With X remaining at unity, deformation 
can occur in the axisymmetric mode; aq.d since = 0 , can take on either 

positive or negative values. If as shown in Figure 2 g^ attains the value 
- 9n/^o f fi n’ pr»t- nf in formation ("18^ vanishes and a bifurcation 

Sm 

of the solution occurs. The bifurcated solution corresponds to falling values of 
X with deformation in both the g^ and g 2 modes. If the deformation is such 
that attains the value 2p/3c , then the coefficient of I3 equation (19) 

vanishes and X falls with deformation in the g^ and g3 modes. In either 
case, the maximum value of X attained is the classical value, X = 1 . 

AXISYMMETRIC IMPERFECTION 

If the imperfection is purely axisymmetric (g 2 =0, g^ > 0) the prebuckling 
deformation is also purely axisymmetric. From (17) 

h * • rh: 
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and P =F =0 until there is a bifurcation of the solution which occurs when the 
2 3 

coefficient of i- n (18) 

1 + P - ^ + f c| L - cl x 

vanishes. Following bifurcation the value of X falls with deformation occuring 
in the ^ as we ll as the axisymmetric mode. This behavior is also depicted in 
Figure 2. The maximum value of X , denoted by X ^ , occurs at bifurcation 
and is found from the previous two equations to be 

2(1 + p -X. M )(1 "V ' C 2 +X M )cI 1 = 0 (20) 

The buckling load for an unpressurized shell is given by 

2(1 - \°) 2 - (2 +k°)c| 1 = 0 . (21) 

and is shown as the dashed curve in, Figure 3 where it can be compared with an upper - 

bound solution obtained in a later section. . vs p curves obtained from 

M 

A •- ~ ~ uoe n-f t hppn chosen such 

^ s O— SJ> - - ^ 

that the unpressurized shells buckle at X - .3, .5 and .7 • These curves as 
well as all others in this paper are calculated with V = 1/3 . 

^ASYMMETRIC AND AXISYMMETRIC IMPERFECTIONS ^ 

If both and are non-zero, deformation occurs in both the and 

§2 modes for any non-zero values of p and X • Non-zero values of can 


44 * 

If < 0 , the coefficient of in (19) vanishes when X reaches X M 

and the bifurcated solution corresponds to falling values of X with 
deformation in the and modes. X ^ is thus dependent only on the 

absolute value of * 
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occur only if the coefficient of in equation (19) vanishes, which for the 

moment we will assume does not happen. 

An expression relating \ and | 2 can be obtained by eliminating from 

the two equilibrium equations (17) and (18). This is 


\ 2 (l 2 + I 2 > " ^ 2 (2 + P> + I 2 (1 + P> + 2 c ^ 2 + C I 2 )1 


+ (1 + p)£ 9 - ( | cg 2 + 4 2 )( cg 2 + | c| 2 g 2 ) --eg,!, + pg, - 0 


>2 b l 
( 22 ) 


A maximum value of X (if it occurs) is associated with = 0 • 

This condition along with (22) yields an expression for X analogous to (20) 


but considerably more complicated. It was found more convenient to obtain X 


M 


from a plot of 1 vs, using (22) . If g^ is positive and large compared 

to I2 j £>1 assumes only negative values and the coefficient of g^ in equation 
(19) 


(1 + P - X) - f c| x + eg 


(23) 


does not vanish. However, when g^ is zero or small compared to negative 

values of g^ are possible; and for sufficiently large values of p it was 
found that this coefficient vanished for a value of X less than the 

M 

provided by (22). Furthermore, in these cases it was found that the bifurcated 
solution resulted in X decreasing with deformation in the g^ as well as 
g^ and g^ modes. In such cases the maximum value of X attained is that 
for which (23) vanishes. 

The indicated calculations were carried out for several combinations of 
g^ and ; and the maximum value of X, X ^ , was obtained over a range 

of p from 0 to 4. X vs p curves are shown in the three plots of Figure 5. 
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In the first plot the combinations of and are such that the unpressurized 

cylinders buckle at , while the other two plots depict cylinders which 

buckle at .5 and when free of internal pressure. The curves of Figure 4 are 
also included. 

It is clear that prediction of the buckling load of a pressurized cylinder 
requires knowledge of the relative amount of axi symmetric imperfection of the 
wave length R/2jtq Q . In effect, the asymmetric imperfections are ironed out by 
the pressure while the axisymmetric ones are not. If i-^/^ i- s sma ll is 

almost the classical value when p is near unity. If, however, the initial 
imperfection is purely axisymmetric the buckling load is much less influenced by 
internal pressure as indicated by the curves = 00 • 

AN UPPER BOUND FOR THE CASE OF AXISYMMETRIC IMPERFECTION 

The previous approximate analyses bares the role of the axisymmetric 
imperfection of wave length R/2itq^ . Not only does it have the largest degrading 
effect for the unpressurized shell, for buckling with p > 1 it almost completely 
determines the buckling load. With relatively little difficulty an upperbound 
to the buckling load as predicted by the non-linear Donnell-shell equations can 
be obtained for a cylindrical shell with axisymmetric imperfections. This is of 
particular interest in light of the role of such imperfections and the approximate 
nature of the previous calculations. 

Only the steps of this calculation will be recorded here since the details 
of a similar calculation for unpressurized shells have been t given by Koiter (8) . 

If the imperfection is of the form 

W = - f-h cos(q x/R) , 
o b l Nn o ’ 
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the state of deformation of the shell can be written as 


W = < + (1 + v ) - T-^-T cos(q x/R)l + w 

l q o ' J 


F = < 


,^Z + i pRx 2 + Y g i h - 

q Q q o (l-M 


cosCq^x/R ) } + f 


(24) 


(25) 


where the terms in the brackets are the exact prebuckling solution to the non- 
linear Donnell -shell equations (1) and (2). Prior to buckling, w and f are 
zero. Buckling occurs with the first deviation from the axisymmetric prebuckling 
deformation. Thus, we look for the value of the load parameter X at which the 
non-linear shell equations admit non-zero solutions for w and f . Substituting 
(24) and (25) in (1) and (2) and linearizing with respect to w and f (we are 
looking for infinitesimal deviations from zero) we find 



w + 
,xx 


R 2 (l - X) 


cos(q x/R) w = 0 

° ,yy 
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cos(q x/R) f 
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+ vT cos(q x/R)w 

yy R(1 -X) . 


yy 


+ 2E^ w - P Rw = 0 (27) 

2 .xx r .yy 


An exact solution to these equations is not found; however, an approximate 
method is employed which guarantees that the estimate of the buckling stress is 
an upperbound to the exact buckling stress. Assume 


w = £h cos(q o x/2R)cos( Yq Q y/2R) 


(28) 
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Solve the compatibility equation (26) exactly for f in terms of the assumed 
w . f must be the sum of terms periodic in the axial and circumferential 
direction since we have ignored the end conditions and the average stresses in 
the shell are represented by the polynomial terms in F in equation (25). Then 
apply the Galkerin procedure to equation (27). This method of solution is 
equivalent to an approximate minimization of the second variation of the energy 
by the Raleigh-Ritz method, which ensures that the eigen-value so obtained is an 
upperbound to the buckling stress.* This straightforward calculation yields the 
following eigen-value equation 

(1 - M 2 [ (1 4 1 * + + 2 P* 2 " -*•)(*.+ yn 

(1 + Y ) 1 1(1 + Y ) 



+ Wip 2 ! 


(1 


. 2,2 

+ Y ) 


(9 


. 2,2 

+ Y ) 


] = 0 


(29) 


The approximate buckling load, , for any given values of p and 

corresponds to the value of y such that X given by (29) is a minimum. 

With p = 0 , (29) is the equation Koiter used to plot an upperbound to the 

buckling load of an unpressurized shell with axisymmetric imperfections. This 
curve for v = 3 is shown in Figure 3. For small values of the minimum 

value of \ from ( 29 ) corresponds to values of y in the neighborhood of unity 
and (29) becomes 

2(1 - X .") 3 - (1 - V <2 + + II ( C K )2 = 0 

An elaboration of this statement is found in Appendix I. 


( 30 ) 
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Equation (2f) from the more approximate analysis is obtained if the term 

in (30) is neglected; and thus as seen in Figure 3, the two curves approach each 
other for sufficiently small . 

For non-zero values of p the minimum X with respect to y is found 
from (29) with the use of a digital computer. Three typical buckling load - 
pressure curves are shown in Figure 6. as obtained from (29) or this figure 

is an upperbound to the buckling load for a shell with the initial imperfection 

and internal pressure corresponding to the parameter p . The three curves 
of Figure 4 obtained on the basis of the approximate calculation do not differ 
appreciably from the upperbound curves of Figure 6, although the value of f 
associated with each approximate curve is smaller than the value associated with 
the corresponding upperbound curve. Use of the curves of Figure 5 would seem 
justified if one wishes to relate the buckling load of a pressurized cylinder to 
the buckling load of an unpressurized cylinder with similar initial imperfections. 

The elastic buckling description is valid only as long as no plastic 
deformation occurs prior to buckling. The elastic hoop strain due to internal 
pressure alone is of the order ph/R . Thus, in most cases plastic deformation 
will occur well before the pressure parameter attains the value 10 ♦ Clearly, 
the buckling load can remain below the classical value for the entire range of 
elastic buckling. 

COMPARISON WITH SOME EXPERIMENTAL RESULTS 

Recent tests with Mylar cylinders performed by Weingarten, Morgan and Seide 
(9) are particularly suited to comparison with the results obtained here. Mylar 
is able to suffer reasonably large strains prior to deforming plastically. Thus, 
Weingarten, et. al., were able to perform a series of tests at different internal 
pressures with the same specimen without incurring any plastic deformation. 
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Figure 7 presents two typical tests series and two theoretical curves chosen 
from Figure 5 which best fit the experimental data, Weingarten, et. al., did not 
report any information with respect to either the form or magnitude of the imper- 
fection which would permit us to assign values to g ^ and • Certainly the 

imperfection representation (12) assumed in the analysis could represent the true 
imperfection only in an average sense; and especially for p < 1 , a more exact 
description would require additional asymmetric terms. Nevertheless, the form of 
the experimental results can be reproduced by an appropriate a posteriori choice 
of and . 

The radius -thickness ratio of the previously mentioned tests ranged from 
200 to 2000, with the maximum load of the unpressurized shells ranging from about 
.6 of the classical value at R/h = 200 to .3 at R/h = 2000. This R/h dependence 
is most readily interpreted in light of the present analysis by associating larger 
imperfections (relative to the shell thickness) with larger values of R/h . 

Indeed, it seems reasonable that such would be the case. 

Judging from the results presented, the relative amount of the axisymmetric 
imperfection is small compared to the asymmetric imperfections in the Weingarten, 
et. al., test specimens. This is particularly the case for specimen 100.2, and 
the buckling load is only slightly below the classical value for p > 1 . Included 
in Figure 7 is data from a series of tests on axially loaded, pressurized aluminum 
cylinders performed at the D.V.L. and reported by Thieleman (2). In this series 
of tests, new specimens had to be used for each test; and although the radius and 
thickness were unchanged, there Was undoubtedly some variation in initial 
imperfections from specimen to specimen as indicated by the data scatter • The 
important feature of these tests is that the buckling load remains well below 
the classical value for values of the pressure parameter well above unity. 
Axisymmetric initial imperfections are strongly suspected. 
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APPENDIX I 


We have asserted that X as obtained from equation (29) is an upper bound 

to the exact buckling load as predicted by equations (26) and (27). The following 

demonstration follows from Koiter's (6) general theory of elastic instability. 

The trivial (or prebuckling) deformation at a given value of load X is denoted 

by w and is construed in a generalized sense. If w is any kinematically 
s 

admissible deformation in addition to the prebuckling deformation w , then the 

s 

change in potential energy of the structure is 

P[w + w g ] - P [w g ] = P 2 [w,X] 4- P 3 [w,X] + 

where P^fw^X] is a functional which is homogeneous and of nth degree in w . 
P-|[w,X] is zero since w is an equilibrium configuration. We have assumed 
w is sufficiently small to make this expansion meaningful. 

At any value of the load parameter, X , the structure is stable if 

■D 1 n -Pov oil arlmi coi'KI o t.t 

£ 

Thus, if for any \ and admissible w , 

a a 

P 2 ( VV = ° (A_1) 

then obviously X > X where X is the lowest value of X such that there 
a — c c 

exists an admissible w for which P^(w, X) = 0 . 

In the present case the terms in the brackets of equations (24) and (25) 
constitute the exact trivial solution w g . Equations (26) and (27) are equivalent 
to the Euler equations for minimizing P^(w,X) . Since we have solved the 
compatibility equation exactly in terms of the assumed radial component of the 
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shell displacement, the stresses are derived from an admissible displacement 
field. The Galerkin solution to the equilibrium equation is equivalent to an 
approximate Raleigh-Ritz minimalization of P^Cw,^) . That is, if the additional 
admissible displacement which we have assumed is denoted by £w , then the 

SL 

approximate eigen-value X is found by 

a 

ti * 2 « 2 v x .> - 2 *wv - ° 

Thus by (A-l) \ > \ 

J a — c 
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FIGURE 1 FIGURE 2 



FIG. 3. EFFECT OF AXISYMMETRIC IMPERFECTION ON 
BUCKLING OF UNPRESSURIZED SHELL. 









FIG. 5c EFFECTS OF AXISYMMETRIC (£,)AND ASYMMETRIC (f z ) 
IMPERFECTIONS. 




FIG. 6 UPPER- BOUND TO BUCKLING LOAD FOR PRESSURIZED 
CYLINDER WITH AXISYMMETRIC IMPERFECTION. 





